Asymmetric Properties of Heat Conduction in a One-Dimensional Prenkel-Kontorova 
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In this Letter, we siiow numerically that the rectifying eflFect of heat flux in a one-dimensional 
two-segment Frenkel-Kontorova chain demonstrated in recent literature is merely available under 
the limit of the weak coupling between the two constituent segments. Surprisingly, the rectifying 
effect will be reversed when the properties of the interface and the system size change. The two 
types of asymmetric heat conduction are dominated by different mechanisms, which are all induced 
by the nonlinearity. We further discuss the possibility of the experimental realization of thermal 
diode or rectifier devices. 
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Understanding of the phonon transport behavior un- 
der nonequihbrium stationary states is much less than the 
electron in charge conduction. It is a key distinctness be- 
tween them that the former is related to nonequihbrium 
statistic mechanics, while the latter is basically the trans- 
port behavior under thermal equilibrium states. Even 
for the phonon thermal conduction in the simplest solid 
model, such as the Fermi-Pasta-Ulam f3 model, in low 
dimensional cases, it is still an open question to get full 
microscopic dynamics description of macroscopic thermal 
transport behavior, i.e. the Fourier law [l|, 0,01 13 ■ On the 
other hand, the band theory predicts the unidirectional 
transport of the electrons in heterogeneity semiconductor 
junctions, which have very useful applications in modern 
electronics. This property is basically attributed to the 
facts that the electrons obey Fermi statistics and only 
the valence electrons take part in the transport process. 
In contrast, the phonons obey Bose-Einstein statistics 
and all take part in thermal transport. In this respect, a 
phonon, viewed as a collective exciton in crystal lattices, 
is not expected to exhibit intrinsic asymmetric transport 
properties, unless the heat conduction of material dra- 
matically changes when the heat baths are exchanged. 

Recently, however, the asymmetric heat conduction 
in one-dimensional inhomogeneous chains with nonlin- 
ear on-site potentials has been reported via computer 
simulations (5j, (61 i7l,i^]. They proposed some mechanisms 
which allow the heat flux in one direction while the sys- 
tems act almost like an insulator when the heat baths are 
exchanged. These mechanisms are based on a common 
idea that phonon bands of the different segments of the 
chain change from overlap to separation when the heat 
baths are exchanged. The phonon band shift controlled 
by temperature is attributed to the different contribu- 
tion of nonlinearity to the effective phonon spectra at 



different temperatures. The effective phonon spectrum 
can be obtained quantitatively by self-consistent phonon 
mcthodPl and qualitatively by linearization treatment 
of Hamiltoniansj^. On the other hand, the numerical 
results in previous work0, Q show that the rectifying 
effect decreases as the system size and the coupling of 
the interface increase. Nevertheless, the picture of over- 
lap/separation of phonon bands, which treats each of the 
different segments of the chain independently, is totally 
independent of the system size and the coupling of the 
interface. Thus, this picture fails to explain numerical 
results on the decreasing of the rectifying effect. Surpris- 
ingly, as shown by our numerical results in this paper, 
the rectifying effect will be reversed when these parame- 
ters further increase. So, more in-depth understanding is 
needed to assess possible applications of this counterin- 
tuitive phenomenon discovered by computer simulations. 

In this Letter, we implement extensive numerical 
studies about a one-dimensional two-segment Frenkel- 
Kontorova (FK) chain that is exactly the same as the 
model studied by Ref.(6i]. We find that the rectifying ef- 
fect demonstrated in Ref.|^ will be reversed when the 
properties of the interface and the system size change. 
A simple physical picture was proposed to understand 
the reversal of the rectifying effect, and a new but very 
intuitive mechanism for the new type of asymmetric 
heat transport was addressed. We further argue that 
the experimental realization is a hard task, because the 
asymmetric heat conduction dramatically depends on the 
properties of the interface and the system size. 

We consider a one-dimensional chain consisting of two 
segments of N/2 particles. The two segments are con- 
nected by a linear spring of constant kint- The Hamilto- 
nian of the whole system is 



H ^ Ha + Hb + -jkint{xN/2+i - xn/2 - a) 



(1) 
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where Ha and Hb are the Hamiltonian of the left seg- 
ment (A segment) and the right one (B segment), respec- 
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FIG, 1: Total heat current ,J± versus interface harmonic cou- 
pling strength hnt for (a) iV = 100; (b) = 1000; (c) 
N = 2000, 



lively, and both are the FK chains described by 

2 ^ y 

^ = + 2^'y^^+^ - ^» - - 7^:32 27ra;,. (2) 

i—l ^ ' 

We set the mass of the particles and the lattice constant 
m = a = 1, and fix Va = 5, Vb = 1, fcyi = 1 and ks = 
0,2, The main adjustable parameters in this paper are 
hint and N. 

In our numerical simulations we use fixed boundary 
conditions and the chain is coupled, at the two ends, 
with heat baths at temperatures T_|_ and r_ respec- 
tively. We take T+ = 0,105 and r_ = 0,035 through 
our paper. We use Nose-Hoover heat baths and integrate 
the differential equations of motion by using the fourth- 
order Runge-Kutta algorithm as described in ^J, We 
compute the temperature profile inside the system, i,e,, 
the local temperature at site i defined as Ti = m{xi ), 
where () stands for temporal average, the local heat flux 
ji — k{xi{xi — Xi-i))^ and J — Nj the total heat flux. 
The simulations are performed long enough to allow the 
system to reach a steady state where the local heat flux is 
constant along the chain. We denote the absolute value 
of total heat flux from left to right as J+ when the bath 
at higher temperature r_|_ is at the left end of the chain 
and as J_ when the baths are exchanged. 

Figure 1 shows kmt dependence of J+ and J_ in dif- 
ferent system size = 100, 1000, 2000, As k,nt in- 
creases there exist a crossover from the situation with 
J+ > to that with J+ < J_, that is, the rectifying 



effect is reversed. In the situation J_|_ > J_, our nu- 
merical results are in agreement with those reported in 
Ref. Q at the same parameters and the prediction by the 
effective-phonon analysis (as shown by the Fig. 4 in the 
Ref. [5). When the heat flux flows from A to B there 
is a big overlap between the phonon bands of A and B, 
while there is almost a separation between them when 
the baths are exchanged. Only phonons with the fre- 
quencies within the overlap range have contributions to 
heat flux. So J+ > J_ and J_ is near zero. Moreover, 
the overlapping bandwidth between A and B is almost 
unchanged when the two segments are very weakly cou- 
pled. So in this case the J± are only determined by the 
transmission of the phonons crossing the interface, i.e., 
kint- One can expect that J± increase in the same way 
as kint increases. This is confirmed by the numerical re- 
sult that j± ~ kf„t with small k^t and N (as shown 
by the Fig. 2(a) in Ref.0|). In fact, this relationship is 
universal when one consider the transmission of phonons 
crossing the interface between two kind of material that 
are weakly joined by a harmonic spring prl IT^ . In the 
situation J+ < J_, the results are strikingly different. 
J± have sharp increases compared to the former case, 
and almost independent of the kmt as shown in Fig. 1(b) 
and Fig. 1(c). These results are disagreement with the 
above prediction based on the effective-phonon analysis, 
and strongly suggest that phonon bands of A and B mix 
up. It is obvious that the band mixing is due to the in- 
teraction among phonons induced by the nonlinearity in 
the chain. In fact, the effective-phonon analysis in Ref@ 
treated each of the two segments of the chain indepen- 
dently, and it could be based on the condition k^t — > 0. 
Consequently, the prediction, i.e. J+ > J^, is well agree- 
ment with the numerical results for very small k^t (a 
typical value in simulations is 0.05). As kint increases 
the two-segment chain behaves as a whole. The phonons 
of the chain in this situation become different from the 
ones of A and B due to the band mixing. Therefore, the 
individual effective-phonon analysis to each of A and B 
fails to explain the new type of asymmetric heat conduc- 
tion. 

There is the difference between curves of the J+ and 
J_ in Fig. 1(a). J+ reaches saturation as kint increases. 
J- reaches its maximum value at the middle kint value 
with {kA + ks)!'^ = 0.6, but decreases as kint further in- 
creases. Specially, J+ is larger than J_ again when kint 
larger than the value kint — kA + ks = 1-2 as shown 
in Fig. 1(a). These phenomena can be understood as fol- 
low. The role of kint changes to an impurity in the two- 
segment chain when kint 3> {kA + kB)/^. The impurity 
just scatters the high-frequency phonons, and does not 
affect the low-frequency ones. We speculate that, for the 
mixing phonon spectrum, the high-frequency phonons 
have main contributions to J_ and the low-frequency 
ones to J+. As a result, J_ drops due to the impu- 
rity scattering, while J+ is almost unchanged as kint 
increases. The effect of the impurity on the heat con- 
duction reduces with N increasing. So the reduction of 
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FIG. 2: Total heat current J± versus the system size A'^ in the 
log-log plot for (a) ki„t = 0.05; (b) ki„t = 0.2; (c) kint = 0.6. 
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J_ decrease as N increase as shown by Fig.l. It is clear 
seen that both J_(. and J_ almost saturate as ki^t in- 
creases as shown in Fig. 1(c). This means that the heat 
conduction behavior of the chain is almost independent 
of kint when the two segments are well coupled. Note 
also that the crossover occurs at kint ~ 0.5, 0.24, 0.15 
for = 100, 1000, 2000. This result strongly implies 
that the crossover will also appears as N increases under 
certain fixed kint- It is confirmed by Fig. 2. 

Figure 2 depicts the N dependence of J+ and J_ 
at different kint- For a typical weak coupling value, 
kint = 0.05, the crossover does not occur, i.e. J+ > J-, 
up to TV = 3000. It is clear, however, seen that the 
increasing of J_ is more faster than as shown in 
Fig. 2(a). This implies that the crossover will occur at a 
larger TV. In Fig. 2(b), the crossover appears at N ^ 1500 
for kint = 0.2. In Fig. 2(c), the two-segment chain ex- 
hibits J+ < J- for the whole regime of N used for simu- 
lations. This result can also be understood by the role of 
kint in the limit case N — s- oo. Strictly speaking, there is 
not a geometric boundary within a one-dimensional sys- 
tem in thermodynamic limit, unless the coupling of the 
two segments is exact zero. For any finite kint, the inter- 
face of two segments behaves like a impurity in the limit 
N ^ oo, and the effect of the impurity on heat conduc- 
tion of the whole system reduces with N increasing, as 
stated above. This suggest that the two-segment chain 
will exhibit the behavior of J+ < J_ in the thermody- 
namical limit. 

Up to now, it is clear seen that there are two types 
of asymmetric heat conduction in the two-segment FK 
model and the transition between them will occur 
by varying the kint and N. When the two segments are 
very weakly coupled, J+ is larger than J_; while J_ is 
larger than J+ when the two segments are well coupled or 
the chain is long enough. Phonon band shift and phonon 
mixing play important roles in the former case and the 
latter case, respectively. It is interesting that these two 
contrary behaviors are all attributed to the nonlinearity. 



FIG. 3: Total heat current J± versus interface anharmonic 
coupling strength Pint for (a) N = 100; (b) N = 1000. 



In order to further verify the role of nonlinearity in 
the asymmetric heat conduction in model we add a 
quartic nonlinear coupling to the interface in model l|T]l. 
as described by 



H ^ Ha + Hb + -fci„f (xjv/2+i - - af 

+ ^f3int{xN/2+l - Xn/2 - a)'^- (3) 

Figure 3 plots the J± versus (3int for N — 100 and 
N = 1000 at a typical weak coupling parameter kint — 
0.05. By increasing the Pint, one can also see the simi- 
lar crossover as shown in Fig.l and Fig. 2. The crossover 
occurs at Pint — 0.06 and Pint — 0.02 for = 100 and 
1000, respectively. These results confirm that the phonon 
mixing due to the nonlinearity leads to the crossover from 
the case J+ > J_ to the one J+ < J_, even if the non- 
linearity is very small. 

A problem immediately arises from the above analysis: 
why is J_ larger than J_|_ when the phonon bands of A 
and B mix up? In the thermodynamical limit, it can be 
expected that the heat conduction of this two-segment 
chain is independent of the properties of interface and 
the system size and determined by the heat transport 
properties of A and B. In Fig. 4, we plot the tempera- 
ture dependence of the heat conductivity k for A and 
B, respectively. As the temperature change from 0.02 to 
0.12, KA decreases, and kb increase. Thus, both A and 
B have the higher heat conductivity when A is contacted 
with r_ = 0.035 and B with T+ = 0.105. When the 
baths are exchanged, both A and B have the lower heat 
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FIG. 4: Heat conductivity ti versus temperature To for seg- 
ment A and segment B, respectively. The temperatures of 
heat baths are r± = ro(l ± 5). To = 0.02, 0.04, 0.06, 0.08, 
0.1, 0.12, and S = 0.1. k = J/{2STo). 



conductivity. It is very intuitive that J_ is larger than 

In summary, we find the reversal of the rectifying effect 
in the one-dimensional two-segment FK model via com- 
puter simulations. We distinguish two different (even 
contrary) effects of the nonlinearity on the heat conduc- 
tion behaviors. For the weak coupling of the two seg- 



ments, the phonon band shift due to the nonlinearity 
leads to J+ > J_. In an other case, the nonlinearity 
causes the mixing of phonon bands and leads to J+ < J_ . 
These two types of asymmetric heat conduction can tran- 
sit to each other by varying the properties of the interface 
and the system size. 

The asymmetric heat conduction discovered by com- 
puter simulations opens new possibilities to design a ther- 
mal rectifier or thermal diode in theoryHlil. However, 
as shown by our numerical results, the function of such 
a thermal device will dramatically depend on the prop- 
erties of the interface. To exactly control the interface 
properties is very difficult in the labTsj, and this diffi- 
culty really leads to the poor reproducibility of the ex- 
perimental measurement of the thermal boundary resis- 
tance. As a result, the function of the thermal devices 
based on the principles stated here will be unpredictable 
due to the unknown interface conditions, such as the 
weak nonlinear atom-atom interaction. Thus, this is a 
hard task to make a thermal rectifier or thermal diode in 
the lab. 
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